Slow Thermalization Between a Lattice and Free Bose Gas 
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Using a 3D spin-dependent optical lattice, we study thermalization and energy exchange between 
two ultracold Bose gases, one of which is bound to the lattice and another that is free from the lattice 
potential. Disruption of inter-species thermalization is revealed through measurements of condensate 
fraction after the lattice is superimposed on the parabolic confining potential. By selectively heating 
the lattice-bound species and measuring the rate of heat transfer to the free state, suppression of 
energy exchange is observed. Comparison with a Fermi's golden rule prediction confirms that this 
effect is caused by a dispersion mismatch that reduces the phase space available for elastic collisions. 
This result has critical implications for methods proposed to cool strongly correlated lattice gases. 
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Ultracold atom gases trapped in optical lattices are 
an ideal system for exploring strongly correlated quan- 
tum matter. Initial work in this area focused on investi- 
gating equilibrium phases of Hubbard models, including 
the realization of the superfluid-to-Mott-insulator phase 
transition for bosons [ll and the metal-to-Mott-insulator 
transition for fermions [H, Q . A new frontier in this field is 
out-of-equilibrium dynamics; recent results include mea- 
surements of quantum quenches 0, [H[ , expansion dynam- 
ics @ , and conduction [tJ . These measurements provide 
insight into phenomena important to material applica- 
tions, such as transport and diffusion. Understanding 
strongly correlated dynamics is also critical to develop- 
ing new techniques for cooling lattice gases. Presently 
attainable temperatures in optical lattices are too high 
to observe many phases of interest, such as antiferro- 
magnetism New methods for cooling are therefore 
required to reach these low-entropy states. 

Here we examine the dynamics most relevant to 
cooling — thermalization and energy exchange — in a 
unique strongly correlated system: a fully three- 
dimensional, species-specific optical lattice [9]. We work 
with a mixture of two bosonic atomic species, one that is 
strongly correlated and bound to the lattice, and another 
weakly interacting state free from the lattice potential. 
This type of lattice has been proposed as a method for 
realizing exotic bosonic and fermionic superfluid states 
a Kondo lattice model UlEIj, and for quantum 
magnetism [3]. The measurements we present here im- 
pact those proposals and suggestions to use this system 
as a platform for cooling and thermometry of strongly 
correlated lattice gases. For cooling, a large number of 
free atoms may be used as a heat reservoir to directly ab- 
sorb entropy from the lattice gas when it is compressed 
[lH or to accept entropy via band decay processes [l6| , 
while a small impurity of free atoms could be employed 
for thermometry [9|. 

There are two methods for realizing species-specific lat- 
tices. In one scheme, different atomic species or isotopes 
are employed, and the relative lattice potential depths 



are adjusted via the lattice laser wavelength [171. Il8j. We 
use the other technique, which is to employ a mixture of 
atoms in different hyper fine (i.e., spin) states trapped in 
a spin-dependent potential [19j. Spin-dependent lattices 
in ID have been used to generate entanglement [20], real- 
ize quantum walks [21[, create atomic impurity disorder 
[22I ]. and to implement a matter- wave probe [23]. We 
have used a 3D spin-dependent lattice to create mixed 
superfluid and Mott insulator phases [9[. 

We use the \F = l,ra F = -1) ("|1)") and 
|F = l,ra F = 0) ( a |0)") states of 87 Rb atoms con- 
fined in a 1064 nm crossed-dipole trap with a geometric 
mean trap frequency (44.5 ± 0.6) Hz; the full details of 
our apparatus can be found in in Ref. [24[. A mixture 
of 75 x 10 3 |1) and 20 x 10 3 |0) atoms is created from 
a Bose-Einstein condensate (BEC) of |1) atoms via 
adiabatic rapid passage (through the |2,0) state) driven 
by a microwave- frequency magnetic field. This mixture 
is miscible in the dipole trap and stabilized against 
spin-changing collisions by an approximately 10 gauss 
magnetic field [25[. The atoms are transferred into a 
spin-dependent cubic lattice by slowing superimposing 
three orthogonal pairs of counter-propagating, linearly 
polarized A = 790 nm laser beams in a lin-±-lin 
configuration (Fig. [I]). At this wavelength, the lattice 
potential depth is proportional to \rrtF\- Thus, the |1) 
atoms are lattice-bound and realize the Bose-Hubbard 
(BH) model, and the |0) atoms do not experience the 
lattice, as shown in Fig.[TJ Using a combination of direct 
imaging and microwave spectroscopy [25], we confirm 
that the two species are spatially overlapped on the 
trap and lattice length scales for all the lattice potential 
depths explored in this work. Using Kapitza-Dirac 
diffraction of a BEC in the |1) state, we calibrate the 
lattice potential depth s, which we express in units of 
the recoil energy E R = 2h 2 7T 2 /m\ 2 (where m is the 
atomic mass). 

By measuring condensate fraction, we observe that 
turning on the lattice disrupts inter-species thermaliza- 
tion. The lattice laser intensity is turned on over 50 ms, 
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FIG. 1. (Color online) Schematic of the spin-dependent lat- 
tice potential. The |1) atoms (blue) are lattice-bound, while 
the |0) atoms (red) do not experience the lattice poten- 
tial. Confinement of both species is provided by a 1064 nm 
crossed-dipole trap (orange lines). The spin-dependent lat- 
tice is formed by three pairs of linearly polarized counter- 
propagating 790 nm laser beams (red lines) with orthogonal 
polarizations (indicated by dashed arrows); full details of the 
spin-dependent lattice are given in Ref. [3]. 



slower than the observed 30 ms inter-species thermaliza- 
tion time in the dipole trap [9]. The lattice modifies 
the dispersion of the lattice-bound atoms and reduces 
their kinetic energy, and therefore, if the lattice turn-on 
is isentropic, the |1) atoms will cool. The temperature of 
the |0) atoms should likewise decrease if the species re- 
main in thermal equilibrium, and thus condensate frac- 
tion Nq/N for that species will increase. We measure 
No/N for each species using time-of- flight imaging as the 
lattice potential depth s is varied. Separate images of 
each component are obtained after turning off the lat- 
tice and dipole trap and 17-18 ms of time-of- flight; the 
components are displaced from each other using a mag- 
netic field gradient. Condensate fraction for the |0) atoms 
is measured using a two-component fit to images, and 
No/N is measured for the |1) atoms using the procedure 
outlined in Ref. [26|. Before turning on the lattice, the 
|0)-component No/N w 0.45, corresponding to tempera- 
ture T w 45 nK. 

As shown in Fig. [2j No/N for the |1) atoms decreases 
for higher s, consistent with quantum depletion induced 
by interactions. In contrast, No/N is insensitive to s 
for the |0) atoms, remaining unchanged from the ini- 
tial value in the dipole trap across s = 4-16. To in- 
terpret measured quantities such as No/N throughout 
this work, we use thermodynamic calculations based on 
site-decoupled mean-field theory (SDMFT) 0, EH and 
the local density approximation for the |1) atoms and the 
semi-ideal model [29[ for the |0) atoms. Temperature in- 
ferred in this manner from No/N is shown in Fig. [2] for 
both species as s is varied; T for the |1) atoms deter- 
mined in this manner has an unknown systematic error, 



since SDMFT does not include all relevant low-energy ex- 
citations (jjcj. The temperature of the |0) atoms remains 
fixed up to s = 16, while the temperature of the |1) atoms 
decreases by approximately an order of magnitude, indi- 
cating that the species are not in equilibrium. The data 
show a significant deviation from the equilibrium temper- 
ature predicted assuming that the lattice turn on is isen- 
tropic (solid line in Fig. [2]). At all lattice depths, the |1) 
(|0)) temperature is lower (higher) than the equilibrium 
temperature, consistent with the lattice removing energy 
from the |0) atoms and suppression of inter-species en- 
ergy exchange. This measured deviation from the pre- 
dicted equilibrium temperature for both species indicates 
that insignificant thermalization occurs on a timescale for 
which equilibration is nearly complete in the dipole trap. 

To quantitatively characterize inter-species thermal- 
ization, we selectively heat the |1) atoms and measure 
the rate of heat transfer Q to the |0) component. We 
choose this approach primarily because interpreting the 
heat transfer rate requires knowledge of the inter-species 
temperature difference, which is complicated by the ab- 
sence of a verified, simple method for measuring the tem- 
perature of the |1) atoms at low temperature. In fact, a 
primary motivation for this work is using the |0) com- 
ponent as a thermometer [jj. While the temperature of 
the |0) component can be measured straightforwardly, 
the temperature difference between the two components 
cannot be controlled as the lattice potential depth is var- 
ied without knowledge of T for the |1) atoms. 

We overcome this complication by heating the |1) com- 
ponent to infinite kinetic temperature T with respect to 
the quasimomentum degree of freedom. Although this 
may seem to be an unphysical limit, in a single-band 
lattice system T — >• oo is well defined. It does not cor- 
respond to infinite energy because the kinetic energy is 
bounded by 4Dt, where D is the dimensionality and t is 
the Hubbard tunneling energy. For a three-dimensional 
multi-band system, such as the lattice in our experiment, 
we define the T — >• oo limit as 12t <C /c#T <C E^, where 
Ebg is the bandgap, and ks is Boltzmann's constant. In 
this T — >• oo limit, all quasimomentum states in the low- 
est energy band are equally occupied. We use a well- 
established procedure for generating this configuration — 
a sequence of rapid increases in s known as dephasing 
(3lj . To dephase, we transiently raise the lattice depth 
to s = 20 three times. Because the |0) component does 
not experience the lattice potential, the dephasing pro- 
cedure has no direct effect on their temperature. Images 
of the |1) atoms after dephasing show that T oo has 
been achieved and that the temperature associated with 
the density distribution is unaffected (iBj . 

To measure the heat transfer rate, we turn on the lat- 
tice over 50 ms, dephase the |1) component to infinite 
T, wait a time thoid for heat to transfer to the |0) com- 
ponent with the atoms held in the lattice and trap, and 
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FIG. 2. (Color online) (Top) Condensate fraction versus lat- 
tice depth for each species. Blue (red) points are for the |1) 
(|0)) atoms. Each point is averaged over 9-15 experimental 
runs; the error bars represent the standard error in the mean. 
For data shown as circles, No/N was measured from time-of- 
flight images obtained after turning off the lattice potential 
in less than 50 ns. For data shown as triangles, No/N was 
measured after bandmapping the lattice in 100 /is, and there- 
fore these points may overestimate condensate fraction and 
underestimate temperature 26]. (Bottom) Temperature in- 
ferred from the measured No/N and number of atoms; the 
statistical uncertainty T is less than 6%. Mixed Mott in- 
sulator and superfluid phases exist (at zero temperature) in 
the lattice for s > 13. The lines are predictions assuming 
that entropy is preserved during the lattice turn on. For the 
dashed lines, inter-species thermalization is assumed to be ab- 
sent, and entropy is conserved separately for each component. 
The predication shown by the solid line assumes inter-species 
thermal equilibrium. 



then measure T for the |0) component using time-of- flight 
imaging. The rate of change of T for the |0) component 
Q, determined by varying thoid at fixed s, is proportional 
to the rate of inter-species heat exchange. Direct com- 
parison of Q between different lattice potential depths 
is possible since the temperature difference immediately 
after the dephasing step is independent of s. The tem- 
perature for the |0) component after loading the lattice 
is insensitive to s, and we fix T oo for the |1) compo- 
nent. Furthermore, since collisions transfer only kinetic 
energy between species, Q is solely sensitive to the kinetic 
energy degrees of freedom. 



Sample data for s = 4 are shown in Fig. [3l We observe 
that T for the |0) component increases linearly in time 
after the dephasing step for all lattice depths sampled 
here, which is consistent with the short-time limit of ex- 
ponential relaxation to equilibrium. We define Q as the 
slope of a linear fit of T for the |0) component vs. thoid- 
To verify this technique, we measure Q at different lattice 
depths without the dephasing step and without the |1) 
component present. As shown in Fig. [3j Q is consistent 
with zero for these measurements at all s sampled here, 
demonstrating that any heat transfer observed is induced 
by dephasing and arises from inter-species interactions. 

As shown in Fig. 2J Q decreases monotonically with 
increasing lattice depth and is consistent with zero for 
s > 12. Thus, the inter-species thermalization rate is 
suppressed as the lattice depth increases, in agreement 
with the disruption of equilibrium evident during turning 
on the lattice (Fig. [2j). While Q is expected to decrease 
with s since the total energy added to the |1) atoms dur- 
ing dephasing is proportional to t (which approximately 
decreases exponentially with s), energetic arguments can- 
not explain Q disappearing for s > 12. 

We compare the measured Q to a Fermi's golden rule 
(FGR) calculation for the inter-species energy exchange 
rate. We treat both states as uniform gases occupying a 
volume V, with the |1) atoms as non-interacting lattice 
gas at T = 00 and the |0) atoms as a weakly interacting 
Bogoliubov Bose gas at T = 0. Inter-species energy ex- 
change proceeds via collisions in which a lattice particle 
of initial quasimomentum q scatters from the A^-particle 
condensate at rest into q 1 and creates a Bogoliubov ex- 
citation of momentum p. The rate for this process is 

2tt , 



T = 



W,?\VM 0>r S [E(q) - E(q') - c(p)] , (1) 



where E(q) = 2£^L 1 [1 — cos(g^(i/fo)] is the energy of 
the |1) particle in the lattice (i indexes lattice directions), 

and e(p) = y (cp) 2 + (p 2 /2m) 2 is the energy of a Bogoli- 
ubov excitation in the |0) gas (c is the speed of sound 
in the |0) condensate). The interaction Vj nt is the con- 
tact interaction for s-wave collisions between ultracold 
atoms summed over all 7V^ ^ |0) condensate particles. 
The matrix element in Eqn. [T] is \{q' ,p \V\nt\ <f, 0)| 2 = 

7V( ° ) (^F") 2^fn) S M'^ where a is the scattering 



length and the Kroenecker-delta represents momentum 



conservation 
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3|. To calculate the total rate of en- 
ergy transfer E = ^Te(p)p^ to the |0) atoms, 
we sum over all final excitation momenta p and average 
over a uniform distribution of initial quasimomenta, 
where is the number of |1) atoms in the volume V. 
Finally, we apply the local-density approximation and in- 
tegrate over the trapped condensate density profile, and 
we approximate the |1) density as one particle per site. 

In qualitative agreement with the measured Q, the 
predicted E (dashed line in Fig. 4) decreases with in- 




FIG. 3. (Color online) (a) Condensate fraction versus hold time for s = 4. Each point is the average of 10 experimental runs, 
and the error bars show the standard error of the mean, (b) Temperature inferred from the data in (a). The line shows a linear 
fit, the slope of which is Q. (c) Control data for the heating rate. The blue circles show Q when the lattice atoms are not 
dephased, and the red squares show Q when no lattice atoms are present. 
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FIG. 4. (Color online) Energy exchange rate for different lat- 
tice potential depths, all in the superfluid regime of the BH 
model. Each point is derived from a fit to data at constant 
s, such as those shown in Fig. 3a; the error bars are the 
fit uncertainty. The lines show a FGR calculation of the en- 
ergy exchange rate E\ the scale has been adjusted so that the 
rates at s — 4 approximately coincide. Averaging over the 
|1) density profile does not significantly alter these results. 
For the dashed line, energy is conserved in collisions, while 
for the solid line the delta- function in Eqn. [1] is replaced by 
a Gaussian with a standard deviation of U . The inset shows 
the fraction of phase space / for elastic collisions, which has 
been averaged over the |0) density distribution. 



creasing lattice potential depth and vanishes above a 
threshold. The essential ingredient in the FGR calcu- 
lation that gives rise to a threshold lattice depth for 
thermalization is conservation of energy and momen- 
tum in collisions. Since the |1) and |0) atoms have en- 
tirely different dispersion relations E(q) and e(p), sat- 
isfying both energy and momentum conservation can 
reduce the phase-space of allowed collisions that ex- 
change energy [l6[. The inset to Fig. 4 shows the frac- 



tion / = £~ p 9 5 [E(q) - E(q>) - c(p)] 5 M , +f / £-1 of 
phase space (averaged over the |1) quasimomentum dis- 
tribution) for elastic collisions, which vanishes for s > 9. 
The observation of a threshold in Q is strong evidence 
that dispersion mismatch prevents thermalization be- 
tween the |1) and |0) gases. The same mechanism is 
responsible for Kapitza resistance [34], which occurs at 
inter facial surfaces. 

Evident in Fig. 4 is a disagreement between the ob- 
served and predicted thresholds. To test if the strong 
interactions between |1) atoms may be responsible for 
this discrepancy, we heuristically account for interactions 
in the FGR calculation by relaxing energy conservation. 
The solid line in Fig. 4 shows E when the J-function in 
Eqn. [D is replaced by a Gaussian with a standard devi- 
ation equal to the Hubbard interaction energy U. The 
close agreement between this case and the measured en- 
ergy exchange rate suggests that a full understanding of 
inter-species thermalization will require a more sophisti- 
cated treatment of the strong interactions induced by the 
lattice. 

The suppression of energy exchange we observe has 
critical implications for schemes that rely on fast energy 
exchange between a lattice and a free-particle bath. Im- 
purity thermometry ^ and cooling methods that require 
equilibration [l5| are problematic, while techniques that 
rely on thermal isolation may be viable. In particular, 
our measurements support the feasibility of the cooling 
technique proposed in Ref. [16[, which relies on one-way 
heat transfer to a free-state bath. Further measurements 
on dynamics in species-specific optical lattices may also 
provide new insight into open questions regarding inter- 
facial effects, such as Kapitza resistance. 

We acknowledge Erich Mueller and Stefan Baur for 
suggesting the T = oo technique and for initial guidance 
on the FGR calculation. This work was supported by the 
DARPA OLE program. 
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